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ABSTRACT 

Large  single  crystals  of  benzene  (C5H5)  were  aligned  by  an 
optical  method.  The  aligned  crystals  were  pulsed  with  ultrasonic 
waves  at  a  frequency  of  approximately  10  megacycles  per  second  using 
quartz  crystal  transducers.  The  velocity  of  the  ultrasonic  waves 
was  measured  by  comparing  the  echo  transit  time  with  that  of  simul- 
taneously generated  waves  moving  through  a  variable  length  path  in 
a  medium  in  which  the  velocity  of  sound  was  known  precisely. 

Velocities  of  longitudinal  and  transverse  waves  along  each  of 
the  principal  axes  of  the  orthorhombic  crystal  were  used  to  compute 
the  main  diagonal  terms  of  the  C^j  matrix.  The  off -diagonal  terms 
were  then  determined  by  velocity  measurements  along  arbitrarily 
selected  directions. 

Velocity  measurements  were  made  and  the  elastic  constants 
computed  over  the  temperature  range  170°K  to  250°K. 

Absorption  in  the  crystal  was  measured  over  the  same  temperature 
range  for  each  of  the  directions  used  for  velocity  measurements.  The 
absorption  was  determined  to  be  aelotropic. 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  0.  B.  Wilson  of  the  U.  S. 
Naval  Postgraduate  School  in  this  work. 
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INTRODUCTION 

Benzene,  the  fundamental  aromatic  hydrocarbon,  possesses  an 
orthorhombic  crystalline  structure  in  the  solid  state. 

The  low  melting  point  of  benzene  (3>.2°C)  renders  the  use  of 
standard  diffraction  methods  for  the  determination  of  orientation 
difficult.  Benzene  crystals  are,  however,  transparent  to  visible 
light  and  possess  a  high  positive  birefringence;  thus  their  optical 
properties  provide  a  suitable  means  of  alignment  using  simple  optical 
methods . 

The  formulation  of  the  strain  energy  equation  for  an  orthorhombic 
crystal  results  in  nine  independent  elastic  constant s(l) .  By  measure- 
ment of  the  velocities  of  longitudinal  and  transverse  ultrasonic  waves 
along  the  three  principal  axes,  the  six  main  diagonal  terms  of  the  C^i 
matrix  can  be  computed.  The  remaining  three  elastic  constants  must 
then  be  determined  by  velocity  measurements  along  arbitrarily  selected 
directions,  using  the  previously  determined  values  for  the  main 
diagonal  terms. 

Using  the  techniques  for  alignment,  pulsing  and  velocity  measure- 
ment described  in  this  thesis,  the  velocity  data  was  obtained  for  the 
computation  of  the  elastic  constants  over  the  temperature  range  170°K 
to  2£0°K. 

In  obtaining  the  velocity  data,  the  opportunity  was  taken  to 
observe  and  record  the  absorption  of  the  ultrasonic  pulses  at  a  fre- 
quency of  10  megacycles  per  second  as  a  function  of  crystal  orientation 
and  of  temperature. 


CRYSTAL  GROWING  AND  CRYSTALLINE  STRUCTURE  OF   BENZENE 

The  large  single  crystals,  necessary  in  order  to  follow  the  align- 
ment and  cutting  procedures  developed,  were  grown  in  a  melt  of  reagent 
grade  benzene.  The  melt  was  contained  in  a  sealed  container  which  was 
maintained  at  a  temperature  of  approximately  5°C.  After  degassing,  the 
bath  was  left  at  this  temperature  for  approximately  two  weeks  during 
which  time  the  crystals  formed. 

The  conglomerate  mass  of  crystals  was  then  broken  apart  and  large 
single  crystals  of  irregular  shape  were  obtained  and  stored  at  a 
temperature  of  about  -10°C  until  required. 

Crystalline  benzene  is  of  the  orthorhombic  class,  space  group 
Pbca,  the  unit  cell  having  the  dimensions  (2)s 

a  -   7.292  A° 
b  =   9.U71  A0 
c  =  6.7U2  A0  at  -195°C 
Optically  the  crystal  is  transparent  with  indices  of  refraction  (3) 

>  1.6U 
-  =*.  =  0.00k 
IS  -OL    >   0.12 
The  crystals  were  thus  very  suitable  for  examination  by  optical 
methods. 
CRYSTAL  ALIGNMENT 

In  order  to  obtain  the  necessary  data  to  solve  the  strain  energy 
equation  and  thereby  compute  the  elastic  constants,  crystals  were 
required  which  were  accurately  aligned  along  each  of  the  principal 
axes  together  with  at  least  one  sample  which  could  be  aligned  along 
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some  arbitrary  axis* 

The  method  of  alignment  used  consisted  of  examining  the  raw  ir- 
regularly shaped  crystal. s  under  a  conical  beam  of  polarized  light. 

For  this  purpose,  a  conoscope  was  constructed  in  the  form  shown 
schematically  (Fig*  1).  For  actual  examination  the  instrument  was 
located  in  a  refrigerated  compartment  maintained  at  a  temperature  of 
about  -10°C. 

In  the  raw  state,  since  there  seems  to  be  no  apparent  preferred 
direction  of  growth  for  benzene  crystals,  initial  examination  was  made 
with  the  crystal  immersed  in  a  glass  bottomed  tank  containing  a  mixture 
of  benzene  (C^H^)  and  methyl  salicylate  (CgH^Oo) .  This  mixture  was  in 
the  proportions  of  bhO   parts  methyl  salicylate  to  300  parts  benzene  by 
volume.  The  refractive  index  of  this  mixture  is  approximately  1»5>2, 
close  to  that  of  benzene  at  ~10°C.  The  actual  mixture  proportions  were 
found  to  be  not  critical  .  Changes  occur  due  to  evaporation  and  the 
solubility  of  the  crystalline  benzene  in  the  mixture.  Fortunately 
the  mixture  is  close  to  the  solubility  limit  of  benzene  in  methyl 
salicylate,  so  that  loss  of  crystal  by  dissolution  is  small. 

The  tank  was  positioned  in  the  conoscope  so  that  the  bottom  of  the 
fluid  was  at  the  focal  point  of  the  convergent  lens,  with  the  divergent 
beam  spreading  up  through  the  fluid. 

In  this  beam  the  raw  crystal  was  examined  for  the  location  of  the 
optic  axes,  the  bisectrix  of  which  locates  the  b  axis  of  the  crystal 
(Fig.  8). 

Examination  was  carried  out  by  rotating  the  crystal,  held  in  a 
pair  of  tongs,  in  the  beam  to  determine  two  extinction  planes „  Since 
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the  normals  to  these  extinction  planes  form  an  orthogonal  system  with 
one  of  the  optic  axes,  this  axis  was  visible  at  90  degrees  to  both  these 
normals. 

This  alignment  was  retained  by  transferring  the  crystal  to  a  special 
clamp  (Fig.  2),  designed  to  accommodate  a  range  of  shapes  and  sizes  of 
crystals  and  at  the  same  time  permit  a  more  accurate  alignment.  The 
clamp  consisted  of  a  short  length  (about  5cm)  of  10cm  diameter  aluminum 
tubing  into  which  were  mounted  in  gimbal  fashion  two  aluminum  rings.  On 
the  inside  of  the  inner  gimbal  were  secured  four  strips  of  spring  steel 
bent  into  the  center  forming  a  spring  collar  which  retains  the  crystal. 
Spring  pressure  was  just  sufficient  to  hold  the  crystal  firmly  in 
position. 

Mounted  in  this  clamp,  the  crystal  was  returned  to  the  bath  with 
the  base  of  the  clamp  forming  a  reference  face  resting  on  the  glass 
bottom  of  the  bath  and  was  thus  in  a  plane  normal  to  the  optic  axis  of 
the  conoscope. 

With  the  freedom  afforded  by  the  gimbals,  the  crystal  could  be 
rotated  with  respect  to  the  optic  axis  of  the  conoscope  to  bring  the 
bisectrix  into  coincidence  with  this  axis.  Since  the  crystal  was  held 
firmly  in  the  clamp,  the  interference  pattern  (Fig.  3)  could  be  clearly 
observed. 

The  alignment  was  preserved  by  removing  the  clamp  and  crystal  from 
the  bath  and,  using  the  base  of  the  clamp  as  reference,  melting  paral- 
lel faces  on  the  crystal,  which  were  finally  covered  with  glass  slides. 

A  check  on  this  alignment  was  made  by  mounting  the  crystal  on  a 
goniometer  attachment,  which  could  be  fitted  to  the  conoscope. 


....  .,. 
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PATTERN. 


FIG.    3. 


Small  corrections  to  the  alignment  could  be  made  by  melting  new 
faces  on  the  crystal,  the  corrections  being  computed  from  the  goni- 
ometer angles  using  the  knowledge  of  the  refractive  indices  of  benzene. 

Alignment  to  the  other  principal  axes  was  accomplished  by  rotating 
the  goniometer  90  degrees  in  the  plane  of  the  optic  axes  (100),  or  in  a 
plane  normal  to  the  plane  of  the  optic  axes  (001) .  After  rotation,  new 
parallel  faces  are  melted  on  the  crystal  and  glass  slides  secured  to 
preserve  the  particular  orientation. 

Any  other  orientation  could  be  obtained  by  appropriate  rotation 
of  the  goniometer.  In  all  cases,  a  knowledge  of  the  direction  of  the 
optic  plane  was  preserved  by  appropriately  marking  the  glass  covering 
slides. 
M0UNTING_  OF  CRYSTALS 

For  the  purpose  of  ultrasonic  pulsing,  it  was  necessary  to  mount 
the  crystals  in  such  a  manner  as  to  permit  electrical  connection  of  a 
radio  frequency  source  and  a  receiver  to  a  transducer,  which  was  in 
mechanical  contact  with  the  crystal.  At  the  same  time,  it  was 
necessary  to  ensure  that  the  crystal  orientation  was  maintained  and 
the  effects  of  sublimation  eliminated.  The  crystal  holder  (Fig.  h) 
was  therefore  designed  and  manufactured  to  meet  these  requirements. 

To  prepare  a  crystal  for  mounting  in  this  holder,  it  was  trans- 
ferred to  a  glass  plate  and  frozen  in  place  using  the  preserved  face  as 
reference;  the  desired  axis  was  thus  maintained  normal  to  the  glass 
plate.  The  cylindrical  shape  required  to  fit  the  holder  was  then 
obtained  by  melting  the  periphery  of  the  crystal,  using  metal  cylinders 
maintained  at  room  temperature  to  melt  the  crystal.  When  the  desired 
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cylindrical  shape  was  obtained,  the  holder  was  slipped  over  the  sample 
and  frozen  in  place .  During  this  process,  the  sample  and  holder  were 
maintained  at  refrigerator  temperature. 

With  the  crystal  in  the  holder,  the  end  faces  were  finished 
parallel  to  the  holder  faces  by  polishing  on  fine  emery  cloth.  The 
end  faces  were  then  sealed  with  l/ii  mil  clear  Mylar  for  the  free  end 
and  l/k   mil  aluminum  coated  Mylar  for  the  end  receiving  the  trans- 
ducer, using  glycerine  as  a  sealing  agent.  The  end  seals  were  secured 
in  place  using  the  holder  end  plates. 

The  transducer  was  mounted  on  the  aluminum  coated  Mylar  using 
Nonaq  stopcock  grease.  These  bonding  materials  were  found  satisfactory 
for  the  transmission  of  both  shear  waves  and  longitudinal  waves  over 
the  temperature  range  for  which  data  was  obtained. 

The  transducers  employed  were  0.5  inch  diameter  quartz  crystal. 
X  cut  crystals  were  used  for  the  generation  of  longitudinal  waves; 
Y  cut  for  shear  or  transverse  waves. 
PULSING  OF  CRYSTALS 

For  pulsing,  the  crystal  holder  was  mounted  in  the  pulsing 
assembly  (Fig.  5).  This  assembly  provided  a  secure  support  for  the 
holder  and  at  the  same  time  permitted  connection  of  the  transducer  to 
the  radio  frequency  source.  The  ground  connection  was  made  through 
the  body  of  the  assembly  and  the  crystal  holder,  the  other  connection 
being  made  through  the  spring  loaded  electrode.  The  method  of  mounting 
the  holder  allowed  adjustment  of  the  electrode  pressure  on  the  trans- 
ducer; a  light  pressure  was  found  best  for  the  generation  of  longitudi- 
nal waves  and  a  heavier  pressure  for  transverse  waves. 
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Radio  frequency  pulses  were  generated  by  an  Arenberg  pulsed 
oscillator,  Model  PG  650-C,  connected  by  low  impedance  coaxial  cable 
to  the  pulsing  assembly.  The  return  echoes  after  reflection  from  the 
free  end  of  the  crystal  were  amplified  by  an  Arenberg  wideband  ampli- 
fier WA-600  and  displayed  on  a  Tektronix  Type  RMli5)A  oscilloscope.  A 
Tektronix  Type  RM  181  time  mark  generator  was  used  for  control  of  pulse 
repetition  rate  and  for  external  sweep  control  of  the  oscilloscope. 
Figure  6  shows  a  schematic  arrangement  of  the  pulsing  system. 
MEASUREMENT  OF  VELOCITY 

To  obtain  accurate  measurements  of  the  velocity  of  the  ultrasonic 
waves  through  the  crystal,  a  modification  of  the  method  of  velocity 
measurement  described  in  (h)   was  used.  Figure  7  shows  the  construction 
of  the  velocity  comparator  designed  and  built  for  this  modified  method. 
In  the  velocity  comparator,  pulses  generated  by  the  quartz  crystal 
transducer  immersed  in  the  liquid  medium  were  reflected  off  the  rigid 
boundary  and  received  at  the  transducer.  The  path  length  and  thus  the 
transit  time  of  the  pulse  was  made  variable  by  altering  the  position 
of  the  boundary  by  means  of  the  micrometer  screw. 

In  making  a  measurement,  the  pulse  trains  from  both  the  comparator 
and  the  crystal  were  superimposed  on  the  oscilloscope  and  the  path 
length  adjusted  to  bring  corresponding  echoes  into  coincidence.  Pre- 
cise matching  could  be  obtained  by  observing  the  interference  effects 
as  the  echoes  were  brought  into  coincidence,  since  the  comparator  and 
benzene  transducers  were  pulsed  in  parallel. 

The  velocity  in  the  crystal  was  then  computed  from  the  known 
velocity  in  the  liquid  medium  by  the  relations 

12 


o 

CO 

CD 

1 

2 

< 

iZ 

o 

LU 

2 

-J 

LU 
Q. 

CO 

LU 

u 

o 
O 
CO 

>- 

CJ 

CO 

CO 

CO 

o 

CD 

z 

H 

CO 

3 

o 

-J 

3 

en 

mmm 

Q 

q: 

a. 

o 

■z. 

lU 

< 

_ 

CD 

u. 

LU 

3 

O 

Q_ 

2 

^ 

< 

or 

>- 

>- 

o 

CD 

LJ  _, 

h; 

H  ^» 

CO 

-J 

3 
Q. 

LU  /^ 
CO  — 
CO  LU 

O 
O 
-J 
LU 

or  ■: 

<  CD 

<  w 

> 

o  £• 

o 

LU 
CO 

_J 

3 

a. 

O 

< 

_J 

o 

CO 

o 

»>, 

\s 

(T 

DC 

o 

< 

H 

2 

< 

LU 

LU 

,-. 

z 

•2S 

LU 

K 

o 

13 


p.  connection 
Insulotor 


Coated    Mylar,  "Seal  a 
Contact 
7y Quartz  Transducer 

y^ Reflecting    Boundary 

/A 

Ethyl    Alcohol 

a    Water 


Micrometer, 
Travel       0-2.5  cm 


VELOCITY  COMPARATOR 


GENERAL  ARRANGEMENT 


FIG.  7. 


/4 


Vc 


1  -10 


Vl 


where:  10  =  initial  micrometer  reading  corresponding  to  zero 
path  length. 

lc  =  length  of  path  through  crystal , 

1  =  micrometer  reading  at  coincidence. 

V]_  -   sound  velocity  in  the  liquid  medium, 

Vc  =  sound  velocity  in  the  crystal , 

The  liquid  medium  used  was  a  mixture  of  100  parts  distilled  water 
to  22  parts  absolute  ethyl  alcohol  by  volume.  This  mixture  transmits 
sound  waves  at  a  velocity  given  by  Willard  (5)  which  is  of  a  magnitude 
suitable  for  comparison  with  the  velocities  of  benzene  crystals  and  is 
substantially  independent  of  temperature  in  the  region  of  23>°C. 
Temperature  control  of  the  device  was,  therefore,  unnecessary.  The 
temperature,  however,  was  monitored  and  corrections  made  for  the 
small  variations  which  occurred. 

Corrections  to  the  measured  length  of  the  crystal  were  necessary 
to  take  into  account  the  effect  of  thermal  contraction  of  benzene. 
These  corrections  were  made  from  the  data  on  the  thermal  expansion  of 
benzene  given  by  Cox  (2),  and  by  Euchen  and  Lindenberg  (6),  and  applied 
as  a  percentage  decrease  in  the  length  of  the  crystal. 

The  accuracy  of  the  velocity  measurements  was  dependent  on  two 
dominant  factors s 

(a)  Measurement  of  crystal  lengths  this  length  was  measured  to 
.001  cm  and,  after  correction  for  thermal  contraction,  was  considered 
accurate  to  within  .01  cm  for  the  3  cm  path  length. 
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(b)  Pulse  matching  between  the  crystal  and  velocity  comparator; 
the  two  10  mc  pulses  were  matched  to  within  one  half  of  one  wave  length, 
resulting  in  a  possible  error  of  .0$ ju   sec  in  12  xj  sec. 

The  possible  errors  introduced  by  other  factors  were  not  of 
significant  magnitude  in  comparison  to  these  two  factors.  Based  on  this 
estimate,  the  velocity  error  was  determined  to  be  less  than  1%. 

Measured  velocities,  using  the  same  crystal,  were  reproduced  to 
within  1%   over  the  temperature  range  2£0°K  to  170°K. 
TEMPERATURE  CONTROL  AND  MEASUREMENT 

Velocity  measurements  were  made  over  the  temperature  range  2£0°K 
to  170°K  by  cooling  the  crystal  mounted  in  the  pulsing  assembly  using 
liquid  air.  In  order  to  maintain  an  even  cooling  rate  the  sealed 
pulsing  assembly  was  inserted  in  a  silvered  Dewar  containing  liquid 
air.  Contact  with  the  liquid  air  was  made  through  a  heat  leak,  which 
permitted  control  of  the  cooling  rate. 

Temperature  was  measured  using  a  Copper-Constantan  thermocouple 
placed  in  contact  with  the  crystal  holder. 
TEMPERATURE  RANGE 

The  upper  limit  250°K  was  predetermined  by  the  need  to  preserve 
the  crystal  from  sublimation  and  also  ensure  the  maintenance  of  a 
satisfactory  bond  between  the  transducer  and  the  crystal. 

Bond  failure  occurred  in  the  temperature  range  1^0°K  to  170°K, 
thus  determining  the  lower  limit  of  170°K.  The  precise  nature  of  the 
failure  was  not  ascertained  but  it  was  believed  that  differential  con- 
traction of  the  crystal  and  the  Mylar  end  seal  resulted  in  separation  at 
the  glycerine  bond  producing  a  void  between  the  transducer  and  the  crystal 
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DETERMINATION  OF  ELASTIC  CONSTANTS 

For  a  crystal  of  orthorhombic  symmetry,  Group  3,  the  nine  inde- 
pendent elastic  constants  may  be  arrayed  in  accordance  with  Cady  (1) 
as  shown: 
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0 

0 

0 

0 

c 

c66 
Since  it  was  possible  to  align  and  mount  crystals  oriented  along 

each  of  the  principal  axes,  the  main  diagonal  terms  C-^  of  the  matrix 

were  computed  directly  from  longitudinal  and  transverse  wave  velocity 

measurements  along  each  of  these  axes. 

The  crystal  orientation,  measured  velocity  and  associated  elastic 

constant  determined  by  this  method  were  as  given  in  Table  1. 


DIRECTION 
OF 
WAVE  PROPAGATION 


100 


_Table  1. 
DIRECTION 
OF 
WAVE  DISPLACEMENT 


100 
010 
001 


010 


001 


100 
010 
001 


100 
010 
001 


VELOCITY 
MEASURED 


ELASTIC 
CONSTANT 


Vaa  i  Cll  *  P  Vaa< 

Vab  !  C66  -  P  Vab2 

Vac  jC£5  »  p   Vac2 

Vba  !  066^7^^ 

Vbb  !  C22  *  P  Vbb2 

Vbc  |  %k  m  f  vbc2 


V 


Vcb    |  Q^  -  p  Vcb2 
Vcc     C33  3  P  vcc2 


Note:  For  axial  velocities,  the  first  letter  indicates  the  axis  of 

propagation,  second  letter  indicates  axis  of  wave  displacement. 
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To  determine  the  off-diagonal  terms  C^j,  the  transformation 
equations  (App.  I)  were  used.  Crystal  orientations  were  chosen  so  that 
each  off -diagonal  constant  could  be  computed  in  terms  of  a  measured 
velocity  and  certain  of  the  previously  determined  main  diagonal 
constants,  For  the  orthorhombic  crystal  the  desired  form  of  the 
equations  was  obtained  by  considering  a  rotation  in  an  axial  plane  to 
an  angle  of  U5  degrees  from  two  of  the  principal  axes  (Fig.  8),  This 
produced  three  different  orientations  giving  three  independent  sets  of 
equations,  each  of  which  could  be  solved  for  one  of  the  three  off- 
di agonal  terms. 

The  three  independent  sets  of  equations  were  obtained  by  sub- 
stituting the  appropriate  values  for  the  direction  cosines  1,  m,  and 
n  in  the  equations  as  demonstrated  in  Appendix  I. 

Substitution  of  these  particular  values  reduced  two  of  the  Cj_-j 
terms  to  zero.  The  remaining  terms  when  substituted  into  the  secular 
equation  for  irrotational  plane  waves  traveling  in  the  X/  direction 
gave  rise  to  the  simplified  determinant  solutions 

(C±i  -  ^V2  )  (Cj£  -  pV*   )  -  0^2  -  0. 
This  may  be  reduced  tos 

where  the  C  terms  are  defined  in  terms  of  the  elastic  constants  and 
the  direction  cosines  (App.  I).  This  results  in  an  equation  in  terms 
of  the  longitudinal  velocity  and  previously  determined  elastic  constants 
(Cii): 

cij  =  \L2/?v2  -  *(cii  +  cjj  +  2Ckk)]2  -  KCjj  -  ciJL)2|    k  .  ckk 
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OPTIC   AXES 


(a)     Principal.  Axes   and   Principal    Planes 


(b)    45°     Oil      Orientation. 


PLANES  AND  AXES 


OF  BENZENE  (CgHg) 


FIG  .  8 
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The  crystal  orientation,  measured  velocity  and  associated  elastic 
constants  were  as  given  in  Table  2; 

Table  2 


Direction  of         Velocity     q,       q„ ,     q..      q 
Wave  Propagation       Measured      1J      1:L      JJ 


110 

V110 

C12 

cll 

Cop 

c66 

101 

V101 

G13 

Cll 

c33 

c3^ 

Oil 

V011 

c23 

Cop 

c33 

CUU 

All  of  the  data  pertinent  to  the  computations  and  the  actual 
values  of  the  elastic  constants  have  been  presented  in  both  graphical 
and  tabular  form,  page  21-26. 
ERROR  IN  ELASTIC  CONSTANTS 

The  sources  of  error  in  the  determination  of  the  elastic  constants 
from  the  measured  velocities  could  be  divided  into  two  categories  -  the 
error  in  crystal  alignment,  and  the  error  in  computing  the  elastic 
constants  from  the  velocity  data. 

The  error  in  crystal  alignment  depended  not  only  on  the  accuracy  of 
the  optical  alignment,  but  also  on  the  method  of  preserving  this  align- 
ment during  the  mounting  procedure.  Since  for  all  alignments  other  than 
that  along  the  b  axis,  the  alignment  of  the  crystal  could  not  be 
rechecked  after  mounting,  no  correction  could  be  applied. 

Another  source  of  alignment  error  was  introduced  in  the  alignment 
of  the  shear  transducer  to  produce  shear  waves  along  the  axes.  To 
eliminate  this  error,  only  longitudinal  wave  velocities  were  used  to 
determine  the  off -diagonal  elastic  constants. 
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AXIAL     VELOCITIES 


170 


90  210  230 

TEMPERATURE    °K 


250 


LONGITUDINAL     VELOCITIES 


70 


90  210  230 

TEMPERATURE     °K 


250 


Z3 


ELASTIC  CONSTANTS  OF  BENZENE  (C6H6) 
x  10  °  dynes/cm,^ 

Temperature 

°K       Cn   C22   C33  Chk       C55   C66   C12   C23   C13 


250 
2U0 
230 
220 


5.83  6.56  6.1U  1.53  3.78  1.97  3.90  3.52  U.oi 

6.01  6.99  6.39  1.57  3.99  2.1U  h,06  3.57  U.22 

6.2U  7.33  6.66  1.63  U. 17  2.32  h.3h  3.67  it.Ul 

6.56  7.78  6.90  1.68  U.Ul  2.U9  U.UU  3.68  li.l*3 

210            6.82  8.12  7.13  1.72  U.6U  2.63  h.67  3.70  U.53 

200            7.06  8.kh  7.35  1.78  U.85  2.78  U.76  3-76  U.61 

190            7.33  8.69  7.59  1.83  5.08  2.92  U.92  3.89  U.66 

180            7.60  8.95  7.81  1.89  5.31  3. Oh  5.07  3.9U  U.66 

170            7.88  9.26  8.01  1.95  SS3  3  =  18  5.08  3.85  U.80 


,'u 


ELASTIC      CONSTANTS     (C|j) 


70 


l  i  i 

190  210  230 

TEMPERATURE       °K 
Z5 


250 


ELASTIC     CONSTANTS        (Cjj) 


170 


90  210  230 

TEMPERATURE     °K 


250 


For  the  crystals  not  aligned  on  a  major  axis,  a  further  source  of 
error  was  introduced  by  the  anisotropic  thermal  contraction,  which 
produced  a  continual  minor  change  of  alignment  throughout  the  tempera- 
ture range .  Because  of  other  larger  sources  of  error  present  in  the 
determination  of  the  off-diagonal  constants,  no  correction  was  made 
for  this  error o 

The  specific  volume  of  benzene  was  considered  to  be  0.9li25  cc/gm 
for  the  temperature  range  2£0°K  to  170°K,  based  on  data  summarized  in 
Andrew  and  Eades  (8).  Possible  errors  in  this  data  were  considered 
negligible  in  comparison  with  the  other  errors  present. 

Since  it  was  not  possible  to  measure  any  alignment  errors  present, 
an  estimate  of  the  accuracy  was  obtained  by  comparison  of  velocities 
which  determine  the  same  elastic  constant,  For  the  constants  determined 
by  longitudinal  axial  velocities,  the  error  was  estimated  at  less  than 
2%,   based  on  velocities  obtained  from  different  crystal  samples.  For 
the  constants  determined  by  shear  velocities,  the  comparative 
velocities  obtained  also  differed  by  no  greater  than  2%, 

The  values  of  the  diagonal  elastic  constants  (Cj_j_)  determined  from 
axial  velocities  were  required  in  the  solution  of  the  off-diagonal 
elastic  constants  (Cjm);  thus,  any  error  present  in  these  values  was 
unavoidably  included  in  the  off-diagonal  values.  This  resulted  in  a 
possible  overall  error  in  the  off-diagonal  elastic  constants  of  up  to 
10$?. 


Table 

of  Estimated  Error  in  Elastic 

Constants 

Elastic  Constants 

Velocities  Used 

Error 

Cll,  C22,  C33 
ClUp  C55,  C66 
Cl2,  C13,  C23 

Longitudinal,  Axial 
Shear,  Axial 
Longitudinal,  Non-axial 

less  than  2% 

approximately  2% 

up  to  10$ 

27 


M5ASUREMENT  OF  ABSORPTION 

Although  measurement  of  absorption  was  not  the  primary  objective, 

the  method  used  to  pulse  crystals  to  obtain  velocity  data  was  readily 

adaptable  to  providing  absorption  data. 

Absorption  measurements  were,  therefore,  made  in  conjunction  with 
velocity  measurements  by  comparing  the  amplitude  of  successive  pulse 
echoes.  The  comparison  was  made  by  matching  the  leading  echo  to  the 
output  of  the  time  marker  generator  displayed  on  the  second  channel  of 
the  oscilloscope.  The  time  marker  generator  output  was  fed  through  an 
attenuator  which  permitted  an  adjustment  of  amplitude  in  steps  of  one 
decibel. 

In  measuring  the  attenuation  of  successive  echoes,  the  time  marker 
was  matched  to  the  leading  echo  and  the  attenuator  then  adjusted  until 
the  time  marker  amplitude  was  the  same  as  that  of  the  succeeding  echo. 
The  attenuation  between  echoes  was,  therefore,  given  directly  as  the 
attenuator  setting  in  decibels.  These  measurements  were  repeated  for 
as  many  echoes  as  appeared  in  a  given  pulse  train. 

All  absorption  measurements  were  made  at  a  frequency  of  10  mega- 
cycles per  second  using  0,!?0  inch  diameter  quartz  crystal  transducers. 
The  accuracy  of  the  measurement  was  extended  to  within  one  half  of  one 
decibel  by  logarithmic  interpolation  of  the  superimposed  waveforms. 
DETERMINATION  OF  ABSORPTION 

The  average  of  the  attenuator  settings  was  then  converted  to 
absorption  in  nepers  per  centimeter  by  a  calculation  based  on  the 
measured  length  of  the  crystals 

.,         /      /  _\  -  attenuation  (db)    Q  ,0   neoers 
Absorption  (nepers/cm)  =  — , ,„, '.  x  8.68   p 

2  x  1  (cm)  db 
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The  ultrasonic  absorption  computed  from  the  measurements  of 
pulse  attenuation  has  been  presented  graphically  on  page  30  and  in 
tabular  form  on  page  31. 

All  of  this  data  was  based  on  the  attenuation  between  first  and 
second  echoes  only,  thus  involving  a  total  path  length  of  six  centi- 
meters in  the  crystal .  For  this  length  of  path,  the  effects  of 
diffraction  were  small  and  within  the  errors  of  the  actual  measure- 
ments, thus  no  correction  was  made  for  attenuation  due  to  diffraction, 

The  absorption  computed  varied  from  greater  than  0.6  nepers  per 
cm.  to  less  than  0.1  nepers  per  cm.  the  variation  being  a  definite 
function  of  crystal  alignment  in  relation  to  the  direction  of  ultra- 
sonic wave  propagation.  The  absorption  for  all  waves  propagating 
along  the  c  axis  was  significantly  less  than  for  waves  propagating 
perpendicular  to  this  direction. 

Some  variation  of  absorption  with  temperature  was  observed  for 
certain  of  the  crystal  directions.  In  general  this  variation  showed 
a  decrease  in  absorption  with  a  decrease  in  temperature. 
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ULTRASONIC     ABSORPTION 
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ULTRASONIC  ABSORPTION  IN  BENZENE 


Direction 

of 
ropagation 

Temperature 

Range 

°K 

Attenuation 
db/cm. 

Absorption 
nepers 
per  ciTio 

Mode 

100 

250 

-  170 

5.3 

0.61 

Shear  (001) 

010 

250 

230 

-  230 

-  170 

3.3 
U.3 

0.38 
Q.k9 

Shear  (001) 

110 

250 
230 

-  230 

-  170 

6.0 
3.5 

0.70 

o.Uo 

Long. 

101 

2U0 
225 
210 
195 

-  225 

-  210 

-  195 

-  170 

a. 7 

k.2 

3.5 

3.0 

o.5h 
0.U8 

o.Uo 
o.35 

Long, 

010 

250 
230 
210 
200 

-  230 

-  210 

-  200 
«  170 

a. 6 

U.2 

3.7 
3.0 

0.53 
0.U8 

0.1*3 

0.35 

Long. 

100 

250 

2U0 
220 

-  2l|0 

-  220 

-  170 

3.3 
3.0 
2.7 

0,38 
0.35 
0.31 

Long. 

001 

250 
2l;0 

-  2U0 

-  170 

2.0 
1.7 

0.23 
0.19 

Long. 

Oil 

250 
215 
205 

-  215 

-  205 

-  170 

1.7 

1.5 
1.3 

0.20 
0.17 
0.15 

Long. 

001 

250 
210 
190 

-  210 

-  190 

-  170 

2.3 
1.3 
0.7 

0.26 
0.15 
0.08 

Shear  (010) 
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APPMDIX  I 


=  0 


The  method  developed  by  J,  R.  Neighbours  (7)  for  the  determination 
of  elastic  constants  is  based  on  velocity  measurements  of  irrotational 
plane  waves  propagated  along  a  specified  reference  axis  (X  ),  where  the 
reference  axes  are  oriented  so  that  the  Y'  axis  lies  in  the  001  (XY) 
plane  of  the  crystal . 

Using  the  measured  velocity  data,  the  elastic  constants  are 
determined  from  the  solution  of  the  secular  equation: 
(Cii  VV2)    C{6  C-^ 

cl6    (C66  'P1^   C56 
Cl5       c/6  (C^  ~^>V2) 

where  the  primed  constants  are  defined  in  terms  of  the  crystal  elastic 

constants  by? 

Cll  =    ^Cll  *  2^2m2ci2  *  2£2n2C13  +  milc22  +  2m2n2C23  f  n^C33 

+  UmVc^  4  U£2n2C^  +  W2m2C66 
C66  =     'f2m2a-2Cn  -  2X2m2a-2C12  +     I  2m2a~2C22  +  .-£  2n2a"2C^ 

f  m2n2a"2C^  +   {(2  -m2)2a~2C66 
C55  =    i>n2a"2Cn  -  2,(  2m2n2a~2C12  -  2lVc13  +  mAn2a"2C22  -  2m2n2C23 

+  n2a2C33  +  m2a-2(a2-n2)2C^  +  iV2(a2-n2)2C^  +  i|/2m2nV2C66 
0x5  =  -^na-lCll  -  2£2m2na"1C12  -    C2na"1(a2-n2)C13  -  mW^C^ 

■f  m2na-1(a2~n2)C23  4  n^aC^  +  2m2na~1(a2-n2)C[lj4  +  2  ^2na"1(a2-n2)C^ 

-  h  C2m2na~1C66 
cl6  =  -.^ma-1^!  +    Xma"1!  C2=m2)C12  „  I  mn2  aT1^  +  X  m3a-1C22 

•f    l  mn2a"1C23  f  2  :  mn2a=*lcUIi  "  ^  mn2a~lc55  f  2  ( '-marH -C  2-m2)C^ 
C^  =     6^mna~2C]_x  -    /mna""2(  i'2-m2)C12  -    £mnC13  -  ^m3na"2C22 

i     (mnC23  f     /mna-2(a2-n2)Cmi  -     /'mna-2(a2-n2)C^  -2/mna-20  2-m2)Cg5 

where  a2  *    C  2  +  m2 
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and  the  -(,  m,  and  n  terms  are  the  direction  cosines  relating  the  X 
direction  to  the  crystallographic  axes  (a,b,c). 

For  a  wave  propagating  along  a  crystallographic  axis,  the  direction 
cosines  becorae  one  or  zero,  which  results  ins 


'ij 


0 


cii 

"  Cii 

so  that 

the   soluti 

Dn 

is  in 

terms 

of 

^>V2  ■  Cii. 

Example ! 

for 

/  = 

1 

°U  = 

Cn 

m 

/^2aa 

m  = 

n 

a    0 

C££  = 

c<< 

m 

fV2ah 

c66  s  C66  =  ^V2ac 
For  a  wave  propagating  in  an  axial  plane  at  hS   degrees  to  two 
crystallographic  axes,  the  direction  cosines  become  zero  or  4^/2, 
For  the  110  case,  where  the  wave  is  propagating  at  U5  degrees  to  the 
a  and  b  axes,  and  perpendicular  to  the  c  axis? 

im  m  -  Ji/2  C-li  =  J(Cn  +   C22  +   UC66  +  2C12) 

c66  B  KCn  4  C22  -2C12) 

cl6  -  i(c22  -cn) 


n  =  0 
a  =  1 


C13   s  C£6  s  ° 


thus  the  determinant  becomes: 


ci6 


0 

C16  <c66  "  ^v2)  0 

0  0  (C^  -    /?V2) 

2\(n/.        „w2n        p.  ^2] 


=  0 


or   (C55  -^)     L(Cn  -^VO(C66  -^PVO    -  C16^J       -  0 
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which  has  solutions : 
r?  - 


(a)     £>ld  -     C$$     =  UChh   i  C55) 

(b)   p^  *  |(cn  +  c6'6)  i    [i(cu  -  c6g)2  +  c^2] 

The  solution  for  the  longitudinal  velocity  (V^  )  is: 


^v  (  2  =  i(cu  +  c66)  *  «Cu  -  c/6)2  ^  c^2] 


2  ,  1, 


=  i(Cn  +  C22  f  2C66)  +  J  [(C12  +  C66)d   +  i(C22  -  Cn) 
Separating  the  off-diagonal  elastic  constant 


/ 


2       1, 


c12  -  {[  ■  i(cll  *  C22  +  2C66)]      ~  ^C22  "  Cll)  J 


c66 


results  in  a  solution  for  the  off -diagonal  constant  in  terms  of  the 
longitudinal  velocity  and  the  previously  determined  diagonal  constants. 

The  solutions  for  the  other  off -diagonal  constants  are  analagous: 
for  the  101  Case  C13  =  U^pV/   -   |(Cn  +  C33  +  2C^)J-  i(C33-C11)2|  2  -  C^ 

23  =  (L2^2  ~  ^C22  +  c33  +  2C^)]~  t(C33-C22)2)I  _  c^ 


for  the  Oil  Case  C 


It  may  be  noted  that  since  the  solution  is  of  the  form: 

(Cj^  -f-  Ckk)  =  a  known  positive  constant  (k2), 
this  allows  for  the  possibility  of  two  solutions; 

Cij  =  k  -  Ckk 
or    Cij  =  ~(k  +  Ckk) 
The  choice  between  solutions  was  made  on  the  basis  of  the  stability 
requirements  for  a  crystal  of  o r thorn ombic  symmetry.  From  Neighbours  (9), 
this  requires  that  the  principal  minors  of  the  matrix  must  all  be 
positive.  The  matrix  was  therefore  expanded  into  terms  of  the  principal 


3U 


minors,  and  the  following  four  inequalities  obtained? 

cllc22  >  c12 

cllc33  >  C132 

c22c33  >  c23 

cll(c22c33~C23  )  +  2C12C23C13    c122c33  ¥   C13  C22 
These  inequalities  were  then  used  to  show  that  the  solution  giving  a 

large  negative  value  was  not  acceptable. 
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